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Abstract 

A detailed study of the Rayleigh-Benard convection in two-dimensions with free-slip boundaries 
is presented. Pseudo-spectral method has been used to numerically solve the system for Rayleigh 
number up to 3.3 x 10^. The system exhibits various convective states: stationary, oscillatory, 
chaotic and soft-turbulent. The 'travelling rolls' instability is observed in the chaotic regime. 
Scaling of Nusselt number shows an exponent close to 0.33. Studies on energy spectrum and flux 
show an inverse cascade of kinetic energy and a forward cascade of entropy. This is consistent with 
the shell-to-shell energy transfer in wave number space. The shell-to-shell energy transfer study 
also indicates a local energy transfer from one shell to the other. 

PACS numbers: 47.27.ek, 47.20.Bp, 47.52. +j 
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I. INTRODUCTION 



The Rayleigh-Benard convection (RBC) is an idealized version of the thermal convection 
found in nature. The study is useful for understanding the convective flow in atmosphere, 
astrophysics, crystal growth, etc. in addition to its usefulness in investigating the heat 
transport, pattern-forming instability, chaos and turbulence. In the classical RBC, a thin 
horizontal layer of incompressible fluid confined between two fiat plates is heated from below. 
The flow dynamics in RBC is governed by two dimensionless parameters: Rayleigh number 
R, which is the ratio of the buoyancy and the dissipative forces, and the Prandtl number 
P, which is the ratio of the thermal and the viscous diffusive time scales. The critical 
Rayleigh number R^, which is the value of R at the onset of convective motion in RBC, 
is independent of P. The value of Rc depends on the nature of boundary conditions. For 
thermally conducting boundaries, it is 657.5 for free-slip conditions, and 1708 for no-slip 
conditions on velocity fields. The Prandtl number P affects the secondary and higher order 
instabilities. For non-zero P, the primary instability always leads to the stationary patterns 
of straight rolls 1]. Busse and Clever investigated extensively the nonlinear stability 
of various instabilities in the RBC with no slip boundary conditions. Their main results 
are well summarized in the so called 'Busse balloon' [s]. Experiments by Krishnamurti 4, 
^ showed the transition of two-dimensional (2D) convection to a three-dimensional (3D) 
and subsequent generation of oscillatory, chaotic and turbulent convection. Ciliberto and 
Rubio ^ observed in their experiments travelling waves normal to the roll axis. 

Numerical simulations , llO| have played a significant role in the investigation of 



RBC. The 2D convection problem captures some of the essential features of 3D convection 

nr, 

specially at higher values of P |11LI12|]. In addition, the 2D simulations are computationally 



less expensive compared to the 3D simulations. Curry et al. [9] investigated the 2D RBC 
for relatively smaller values of the reduced Rayleigh number r = R/Rc (60 < r < 290) 
and observed one and two frequency convective flows in addition to the stationary rolls. 
Goldhirsch et al. lOj studied numerically the initiation of convective rolls in closed geometry. 
Scaling [l^ and turbulent behaviour [3] are also studied in 2D. We have simulated the 
convective flow in various regimes and investigated the energy spectra and fluxes. 

We present here the results of direct numerical simulations (DNS) of 2D RBC in a pure 
fluid (P = 6.8) confined between stress free flat boundaries. We have used various resolutions 
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as required by the problem, and investigated systematically convective flow in a wide range 
of reduced Rayleigh number r (1.01 < r < 5 x 10^). We observe ordered states (both 
stationary and oscillatory), 'travelling rolls' moving chaotically in a direction normal to the 
roll axis, recurring ordered states, and the soft-turbulent regime. The Nusselt number scales 
with e = r — 1 as e*^'^'^. We have studied the energy spectra and fluxes for both velocity 
and temperature fluctuations. The kinetic energy shows inverse cascade, while the entropy 
shows forward cascade. This behaviour is also observed in shell-to-shell energy transfer in 
wave number space. We observe a local energy transfer in the shell-to-shell study. 

The outhne of the paper is as follows. In Section II we describe the governing equations 
and the numerical method. Section III contains descriptions of various convective states 
observed in our simulation. Section IV and V contains discussions on energy spectra, fluxes, 
and shell-to-shell energy transfers of velocity and temperature fluctuations. In Section VI 
we discuss scaling of large-scaling modes and Nusselt number with Rayleigh number. The 
last section contains conclusions. 

II. HYDRODYNAMIC SYSTEM AND NUMERICAL METHOD 

We consider a thin extended layer of Boussinesq fluid of thickness rf, kinematic viscosity 
i>, thermal diffusivity k, and thermal expansion coefficient a confined between two stress-free 
and thermally conducting horizontal plates. An adverse temperature gradient /? = AT/d 
is imposed across the fluid layer. The hydrodynamic equations are nondimensionalized by 
choosing the length scale as d, velocity scale as K,/d, and temperature scale as AT — pd. 
The relevant dimensionless hydrodynamic equations for the RBC are given by 



where v = (fi,f2,f3) is the velocity fluctuation, 6 is the perturbations in the temperature 
field from the steady conduction state, R — agPd'^/uK is the Rayleigh number, and P = u/k 
is the Prandtl number, and z is the buoyancy direction. Two-dimensional rolls are assumed 
to be parallel to the y axis. Free-slip and perfectly conducting boundary conditions at the 
horizontal plates imply 



dtv + (v • V)v 

dtO + (v ■ V)^ 



-Vp + RP9z + PV\, 

V3 + v'e, 



(1) 

(2) 
(3) 



V- V 
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V3 = d,,V3 = 9 = 0, at 2; = 0, 1. (4) 

The fields are considered to be periodic along the x direction. 

The system of equations with the boundary conditions are numerically solved using 



pseudo-spectral method [15|. We use Fourier basis functions for representation along the x 
direction, and sin or cos functions for representation along z direction. In this 2D simulation 
the velocity field is considered to be confined in the x — z plane. Therefore the velocity 
component in the y-direction is considered to be zero, i.e. 

vi{x,z,t) = '^^UmOn{t) exp^imkcx) cos{miz), 

m,n 

V2{x,Z,t) = 0, 

V3{x,z,t) = '^^WmOnit) exp{imkcx) sm{nnz), (5) 

where kc = tt/V^. Various grid resolutions, 64 x 64, 128 x 128, 256 x 256, 512 x 512, have 
been used for the simulation. The aspect ratio of our simulation is : 1. Time stepping 
is carried out by the standard fourth-order Runge-Kutta (RK4) scheme with time steps 
ranging from 10~^ to 10"^. For some runs we only excite the large-scale modes by properly 
selecting the initial conditions. The simulation is run till the system reaches statistically 
steady-state. 

We carried out simulation for wide range of reduced Rayleigh number r = R/Rc to explore 
many of the possible convective states for 2D RBC We used free-slip boundary conditions 
on the velocity field. This condition is much convenient to formulate and use. We shall use 
the no-slip conditions on the velocity field in the future. The range of r was taken from 1.01 
to 5 X 10'^ {R = 664 to 3.3 x 10''). We fixed the thermal Prandtl number at P = 6.8 which 
is the typical value for water at room temperature. Note that the convection in low Prandtl 
number fluids is dominated by wavy rolls that is ruled out in 2D simulation. Therefore, we 
focus on large Prandtl number (P > 1) simulation for our 2D investigation. 

III. SIMULATION RESULTS ON VARIOUS CONVECTIVE STATES 

Our simulation of 2D RBC spanning a wide range of Rayleigh number exhibits a variety 
of convective states: stationary rolls, quasiperiodic rolls, chaotic flow, travelling rolls, and 
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soft-turbulent flow. Fig. [T] exhibits the time series for large scale modes for these states for 
various values of r. The results in Fig. [1] have been observed for initial conditions when 
the imaginary parts of all the modes are nonzero, and the real parts of these modes are 
zero. The only exceptions are 6'oon modes, which are always real due to the reality condition. 
When we interchange the real and imaginary parts of the Fourier modes (except for the 6'oon 
modes) in the initial condition, the real and imaginary parts of the Fourier modes in the 
flnal states also get interchanged. 

A. Ordered states 

Figure [1] shows the time series for real and imaginary parts of the complex mode Wioi 
for various values of r. The stationary straight rolls are observed for r < 80. All the Fourier 
modes remain real (imaginary), if we choose them initially to be real (imaginary). The Fig. 
la shows the time series for stationary rolls. The stationary rolls bifurcate to the oscillatory 
rolls for r > 80. All the modes still remain either real (imaginary), if the initial conditions are 
chosen real (imaginary). With further increase in r (r > 125) both real and imaginary parts 
of all the modes develop, even if they are initially chosen as real or imaginary. Fig. ^ shows 
such behavior at r=140. The real part oscillates around the zero mean and the imaginary 
part around a nonzero mean as we started simulation with imaginary values for the modes. 
The time period of the real part of the mode W^ioi is twice the period of the imaginary 
part. For r > 145, both the real and imaginary parts of a mode show oscillations about 
nonzero mean (see Fig. [T]l for r=400 ). This phenomenon continues till r = 660. Beyond this 
Rayleigh number the period of real and imaginary parts develop a new frequency, and the 
system bifurcates to a quasiperiodic state. This is observed for 660 < r < 770 {R = 4.3 x 10^ 
to 5.1 X 10^). Fig. [T^ shows the quasiperiodic time series for r = 700 {R = 4.6 x 10^). 

Fig. [2] shows the projection of the phase space in the $5(iyioi)-Q'(6'ioi) plane. The pro- 
jection of the phase space is given for the same values of r for which time-series is given 
in Fig. [H Notice that time periods of oscillations of 53(wioi) and Q'(^ioi) are the same for 
r = 100, 140. However, the period of ^^(^loi) is half of that for the $5(^101 ) for r = 400. All 
the oscillatory states observed between 80 < r < 660 may be categorized in three types of 
different oscillatory states. Fig. [2^ shows the temporal quasiperiodic behavior at r = 700; 
quasiperiodic behaviour is observed for 660 < r < 770. 
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B. Chaos and travelling rolls 



For 770 < r < 890 {R = 5.1 x 10^ to 5.9 x 10^), we observe that the time series of the 
large scale Fourier modes is chaotic (see Fig. [T]F). The real and imaginary parts of a mode 
show small fluctuations around a mean value for long time; and suddenly they flip their sign. 
Whenever the real (or the imaginary) parts of the modes change sign, the imaginary (real) 
parts also change sign. The imaginary and the real parts can also show large fluctuations 
without a flip. The temperature modes closely follow the behaviour of the velocity modes. 
The projection of the phase space in S5(wioi)-Q'(^ioi) plane is shown in Fig.[2]F. An interesting 
consequence of the flip of sign of real and imaginary parts of the modes is the travelling rolls 
instability [l^ . A change of sign in the real and imaginary parts of the modes is equivalent 
to a change of phase of the complex amplitudes. The relative change in the phase of rolls 
induces a lateral movement of the 2D rolls. The convection rolls travel in a direction normal 
to the roll axis. The duration for the lateral motion is very small (~ 0.1 diffusive time scale) 
compared to the duration for which the rolls do not move laterally. The 'travelling rolls' 
instability is shown in Fig. [31 The three frames are taken just before the rolls start moving 
{t = 1.4), during the motion (t = 1.43) and just after the motion(t = 1.5). By the time a roll 
travels by half of its wavelength (vr/fcc), the flow direction reverses globally in the box |l^. 



C. Recurring ordered states 

As we increase r beyond r = 900, the chaotic states cease to exist; instead we observe 
a time independent state. This stationary state is seen to persist for a wide range of r 
(900 < r < 5000). However, this state is qualitatively different from the time independent 
state observed for r < 80. The spectra for both the kinetic energy and the entropy show 
power law behaviour which is discussed in the next section. For 5000 < r < 10000, we 
observe periodic state. The amplitude of oscillations are very small compared to the mean. 
The frequency of oscillations of the modes for this regime is very large (~ 4 x 10^). Again 
the energy spectra Ev{k) and Eg{k) show power law. The energy spectra show power law for 
all values of r > 770 irrespective of the temporal behaviour of the kinetic energy or entropy. 
Recurring ordered states (see Figs. [If, g and [2]F, g) involve energy distribution over wide 
range of wave numbers with power law behavior. 
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D. Soft-turbulent state 



The soft-turbulent convection appears as r is increased to 5 x 10^(i? = 3.3 x 10^) (Fig. [1]). 
The structure of the flows in the soft turbulence regime agrees qualitatively with the 2D 



simulation for high Rayleigh number ll|. Both the real and imaginary parts of the large 
scale velocity mode Wioi show large fluctuations around the zero mean value. Fluctuations 
in ^{WiQi) and Q'(W^ioi) are of the order of 200 dimensionless units, which is much larger 
than the fluctuations observed in the VTioi mode in the periodic regime. The probability 
distribution of temperature fluctuations in the middle of the box show Gaussian behaviour 

n 

as observed by Castaing et al. [17] in the soft-turbulent regime. We observe breaking down 
of the spatial correlations in the flow which culminates in the turbulent motion. 

IV. ENERGY SPECTRUM AND FLUXES 

The kinetic energy spectrum E^{k) is defined as the sum of the kinetic energy of the 
Fourier modes contained in the wavenumber shell [k, k + 1), i.e. 

EM= J2 (6) 

fc<fc'<fc+i 

Similarly we define the entropy as 

Ee{k)= J2 ^l^(k')r. (7) 

k<k'<k+l 

Figure H] shows the kinetic energy and entropy spectra for relatively small Rayleigh num- 
bers r = 70, 140, 400. For this range of r the kinetic energy and entropy spectrum are both 
exponential, i.e., E{k) oc exp{—ak). The exponential nature of these spectra continues up 
to r = 600. The value of a for kinetic energy and entropy lies in the range of 0.4 to 0.8. This 
result of ours including the values of a are in good agreement with Curry et. al. Qj. The 
above result indicates that only a small fraction of Fourier modes are excited appreciably 
for low values of r. 

For r beyond 600, the spectra of kinetic energy and entropy follow powerlaw behaviour 
as shown in Fig. [5l Note that the 2D RBC exhibits quasiperiodic (660 < r < 770), chaotic 
(770 < r < 890), recurring fixed points and periodic states (900 < r < 10"^), and soft 
turbulent (r around 5 x 10^) behaviour. In the quasiperiodic, chaotic or periodic regime 



the power-law exponents for the kinetic energy and the entropy are around —4 and —2 
respectively. 

After studying energy spectra, we turn our focus on the energy flux or energy cascade 
rate. The energy flux emanating from a wavenumber sphere is defined as the total en ergy 



transferred from all the modes inside the sphere to all the modes outside the sphere 18|. 



We compute kinetic energy flux n„ and entropy flux Ug using the formalism described in 



Verma 19|. The flux is defined as the energy leaving per unit time from the inside of a 
wavenumber sphere to the outside of the sphere. This energy transfer takes place from the 
modes inside (giver) the wavenumber sphere to the modes outside (receiver) the sphere. 



This is computed as, 



k p 



where the truncated variables and are defined as follows 

v>(k) 



(8) 



if |k| < K, 
v(k) if Ikl > 



v(p) if \p\<K, 
v<(p) = <^ (9) 

if IpI > K. 

The p summation in Eq. [8] is the convolution sum. Pseudo-spectral method can efficiently 
compute the flux using the truncated variables v'* and v^. We repeat this process for every 
K for which we need the flux. In a similar fashion we can compute the flux 11^ of a scalar 
fleld 6. Convective turbulence is anisotropic, yet we compute energy flux which is an average 
quantity over all the angles. The computation of anisotropic fluxes is quite complex, and 
they will be computed in future. 

The numerically computed values for various wavenumber spheres are shown in Fig. [6] for 
r = 400, 830, and 10^. We observe negative kinetic energy flux for low wavenumbers that 
indicates the inverse cascade of kinetic energy. The kinetic energy flux for a small band of 
higher wavenumber is positive, but its magnitude is rather small compared to the negative 
flux values. The inverse cascade of kinetic energy at low wavenumbers could be due to the 
two- dimensionality of the flow since 2D fluid turbulence exhibits inverse cascade of kinetic 
energy for wavenumbers lower than the forcing wavenumbers. As shown in Fig. [6] the energy 
cascade of the entropy however is forward for the full range of wavenumbers. 

8 



We also observe that the amount of kinetic energy flux increases with the increase of r. 
We will discuss these issues in Sec. IVII In the next section we will study energy transfer 
from a wavenumber shell to another wavenumber shell. 

V. SHELL-TO-SHELL ENERGY TRANSFER 

Energy flux gives the overall energy lost from a wavenumber sphere due to the nonlinear 
interactions. More detailed picture of energy transfer is captured by another quantity called 
shell-to-shell energy transfer. This quantity is specially useful for quantifying the locality in 
turbulence. We compute the kinetic energy and entropy transferred from shell m to shell n 
(T^^ and T^^ respectively) using the method described in Verma [l^. 

We define the shell{n) as the shell that contains wavenumbers shellradius{n — 1) < 
K < shellradius{n) . The radii of the shells are distributed logarithmically for A < K < 
Maxpossihlejinner jradius /2. Hence the radius of the 4*'' shell to the [N shell — 3)*^ shell 
is given by 

shellradius{n) = 8 x 2'^''-^'^ (10) 
where 8 x 2*(^'*'**^''~^) = Maxpossihlejinner jradius/ 2. 

The wavenumber range of the shells 1, 2, 3, 4, N shells — 1, Nshells are [0,2), 
[2,4), [4,8), [8,8 X 2^*), [Maxpossihlejinner jradius /2,MaxpossihleSnner jradius), 
[Maxpossihlejinner Jradius, oo) respectively. Thus, the effective shell-to-shell energy trans- 
fer rate from the mth f -shell to the nth v -shell (Eq. (23)) can be written as , 

(11) 

Similarly we can compute the shell-to-shell energy transfer T^,^ for the scalar field 9. 

In figure [7] we show the shell-to-shell energy transfer for different regimes of r, namely at 
r = 400, r = 830, r = 1000, and r = 10000. The giver shell-index m is shown along y axis, 
and the receiver shell-index n is shown along x axis. In our colour scheme, red is maximum 
positive energy transfer, while blue is maximum negative energy transfer. We can draw the 
following conclusions from our analysis. 

1. Both kinetic energy and entropy shell-to-shell transfers are local, i.e., the most signif- 
icant energy transfer is to the nearest shell, and the energy transfer to more distant 
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P) 



<p<k 

n 



shells decreases drastically. 

2. For low wavenumber shells, the shell-to-shell kinetic energy transfer from m to m + 1 
is negative, i.e., T^'^i ^ < 0, while T^_^^ = — T^'^^ ^ > 0. It corresponds to inverse 
cascade of kinetic energy for low-wavenumber shells. The sign of energy transfer 
for m > 5 is reversed indicating a forward energy transfer for higher wavenumber 
shells. These shell-to-shell energy transfer results are consistent with the flux analysis 
presented earlier. 

3. For shell-to-shell entropy transfer, we observe that T^^^^ ^ > indicating forward 
entropy transfer. This result is consistent with entropy flux calculations presented 
earlier. 



VI. SCALING OF LARGE SCALE MODES AND NUSSELT NUMBER 

In this Section we will describe the variation of the amplitude of large-scale modes, energy 
flux, and Nusselt number as a function of reduced Rayleigh number. Some of the large-scale 
velocity and temperature modes present in the system are Wioi, 6*101, and 6*002 • In § III] 
we discussed the dynamics of these modes at various reduced Rayleigh numbers. Here we 
describe the evolution of the amplitudes of these modes as a function of e = r — 1. As 
shown in Fig. [HI iW^ioil grows with e as e^^^^. The amplitude of the mode 6*101 has different 

scaling for small e and large e. As shown in Fig. [HI for small e, |^ioi| ~ e^'^"^, but for large e, 

\fl I ,-0.34 
rioil ~ e 

We observe in our numerical simulation that ^002 ~ 0.15 for all r's. A careful investigation 
of the energy equations of the large-scale mode 6*101 nicely yields this value for 6*002- These 
arguments are described below, the energy equation of 6*101 is 

l^ioil' = r^(lOl) + 3? [^toiU-101] - (tt' + kl) l^ioil' , (12) 

where kc = 7r/^/2, and 9ft represents the real part. The terms T^(lOl) represent the nonlinear 
interaction terms contributing to mode 6*101. If we truncate the Fourier expansion of the 
velocity and temperature modes and keep the modes VTioi, 6*101 and 6*002, then the nonlinear 
transfer terms turn out to be T^(lOl) = 27r3ft [6*j[■olWlol^oo2]• 
We compute various terms of Eq. (fT2|) using the numerical data in the steady-state regime. 
We observe that in all regimes (ordered or chaotic) T^(lOl) ~ 27r3ft [6*j['oi'"^ioi6'oo2] within 5- 
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10%. In addition to the above results, we find that the dissipative term (vr^ + k"^) |6'ioi|^ is 
somewhat smaller as compared to the other two terms. If we ignore the dissipative term, 
Eq. ( fT2ll yields 6*002 ~ ^VI^tt) ~ —0.15, which is observed quite prominently in all our 
DNS resuhs. 

In addition to the above analysis of the large-scale modes, we also study the energy flux 
as a function of r. We observe that the absolute value of maximum energy flux (|n^|) in the 
inverse cascade regime increases with r. In Fig. [9] we plot III^,! as a function of r. For r in 
the range of 10'^ to 5 x 10^, |n^| ~ r^'^ to a good approximation. Note that 

TT ^,3 ^3*0.6 ^1.8 /1Q\ 

where ml is the large-scale velocity. The above estimate of |n^,| is in qualitative agreement 
with our numerical result |n„| ~ r^'^. The increase in the inverse energy flux with the 
increase in r possibly strengthens the large-scale structure leading to fixed point or periodic 
solution for the large-scale mode ll^ioi- The maximum value of 111^,1 appears to tapers off 
beyond r > 5 x 10^ where soft turbulence regime begins. This aspect of energy flux is under 
investigation. 

Another global quantity of interest is the Nusselt number which is the ratio of total heat 
flux and the conductive heat flux. For Rayleigh number r < 1, Nu = 1, as only conduction 
takes place in this regime. Far from the threshold, Nusselt number shows a power-law 
dependence on the Rayleigh number. There are various predictions for the exponent some 
of which are described below. The theory by Malkus ^] based on boundary layer stability 



gives an exponent of |. Some experimental results 21j are in agreement with this prediction, 
but other experiments suggest an exponent closer to | [l7|, l22|, |23| . It is generally believed 



that the Nusselt number exponents depend on boundary conditions and Prandtl number. 



Refer to [2J]) for a recent review on Nuselt number scaling. 

In figure [10] we plot the Nusselt number Nu as a function of e = (r — 1). e is the 
measure of how far the system is from the onset. Nu vs. e therefore shows how the Nusselt 
number changes as we increase the temperature difference across the horizontal plates. In 
the oscillatory and chaotic regime we have used the mean value of Nu for this plot. The plot 
shows two distinct scalings for Nu. Close to onset, the Nusselt number increases linearly 
with e. However for e > 1, Nu ~ 7-o.33±o.oi^ observe significant deviation from the fit 
for 100 < e < 600; it is interesting to note that the complex modes become time dependent 
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around r = 125 that lie in this band. 



VII. CONCLUSIONS 

In this paper we have presented a detailed study of the 2D simulation of the RBC 
with free-slip boundaries for reduced Rayleigh number r up to 5 x 10^ {R = 3.3 x 10^). 
The simulation is performed with a pseudo-spectral method. We observe stationary states, 
time-periodic states with different features, quasiperiodic, chaotic travelling waves and soft- 
turbulence for different regimes of reduced Rayleigh number r. We observe that stationary 
and periodic states recur after chaos that indicates complex nature of the strange attractor. 

We investigated the large scale complex Fourier modes VTioi and 6'ioi in these regimes. 
The real and imaginary parts of the complex modes couple with each other through nonlinear 
interaction of the higher order modes. However this coupling is active only when r is raised 
to a sufficiently high value above r = 1. In the chaotic regime, the real and imaginary parts 
of Wioi and 6'ioi show random flips. These flips lead to travelling roll instability. The rolls 
travel normal to the roll axis due to this instability. This may lead to global flow reversal in 
the system. The results obtained are in agreement with the experimental results of Ciliberto 
et al. fl. 

We analyzed energy spectra for all the convective states. For the transition regime (1 < 
r < 600) we observe exponential energy spectra for both velocity and temperature field. This 
result is consistent with that of Goldhirsch et al. [10]. For 600 < r < 5 x 10^ within which 
we observe chaotic state, recurring stationary and periodic states, and soft-turbulence state, 
the energy spectra for both velocity and thermal fields are power law with exponents close to 
-4 and -2 respectively. We have not yet investigated the energy spectra for hard-turbulence 
regime. 

We computed energy fiuxes for various values of r. We observe negative kinetic energy 
fiux for low-wavenumber shells, and positive kinetic energy fiux for intermediate wavenumber 
shells. The inverse cascade of kinetic energy could be due to the two-dimensionality of the 
system. This claim is not conclusive since the kinetic energy spectrum exponent is close 
to -4, not -5/3 as observed for the 2D fluid turbulence in the inverse cascade regime. It 
is possible that the hard-turbulence regime has -5/3 kinetic energy spectral index for low 
wavenumbers (for r > 5 x 10"^). The entropy flux is positive for all r indicating a forward 
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entropy \0{k)\'^ cascade. 

To investigate energy transfers in convective turbulence in detail we also study the shell- 
to-shell energy transfers for both velocity and temperature fields. We observe the shell-to- 
shell transfer to be local for both the fields. For kinetic energy, the transfer is backward 
for smaller wavenumber shells, and forward for the higher wavcnumber shells. However the 
transfer for the entropy is forward. These results are consistent with the flux results. 

We analyzed the variation of large-scale modes VFioi, ^loi, and ^002 as a function of r. 
We find that l^ioi ~ r°-^^, while ^101 ~ r°-^''' for small r and ^101 ~ r~^-^ for large r. We 
observe that ^002 ~ 0.15 for all r. This result has been derived using a numerical input from 
our simulation. We compute the Nusselt number scaling, and find that it scales as r°'^^. 
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Figure captions 



FIG. 1: Time series for tlie real and imaginary parts of the mode Wioi for r = 70 (a), 
r = 100 (b), r = 140 (c), r = 400 (d), r = 700 (e), r = 830 (f), r = 1500 (g), r = 7000 
(li), r = 50000 (i). Prandtl number is fixed at P = 6.8. In (h) we plot 3fJ(wioi)/78.5 and 
$j(wioi)/20 for visual clarity. 

FIG. 2: Phase plots in the Q'(H^ioi) - Q'(^ioi) plane. The values of r and P are same as in 
the corresponding plots of figure [H 

FIG. 3: Travelling rolls at r = 830 {R = 5.5 x 10^) and P = 6.8. Here red repre- 
sents hot fluid, and blue represents cold fluid. The convection rolls travel chaoticlly to 
the left. The three frames are taken at dimensionless times t = 1.4, 1.43 and 1.5 respectively. 

FIG. 4: Kinetic energy spectra -Et>(k) vs. wavenumber are plotted in (a), (c) and (e) for 
r = 70 (64 X 64), 140(128 x 128) and 400 (512 x 512) respectively. In (b), (d) and (f), 
entropy spectra Eg{k.) vs. wavenumber are plotted for the same sequence of r. 

FIG. 5: (a) Plots of kinetic energy spectra E^{]<.) and (b) entropy spectra -^^(k) as a 
function of /c at r = 830. Similar plots for r = 10000 (c and d), and for r = 50000 (e and f). 

FIG. 6: (a) Plots of flux of kinetic energy n„ and (b) entropy Ug as a function of 
wavenumber sphere radii at r = 400. Similar plots for r = 830 (c and d), and for r = 10000 
(e and f). 

FIG. 7: (a) Plots of shell-to-shell energy transfer for v ^ v channel (Z^^) i^) 0-^9 
channel (T^^) a function of m and n at r = 400. In the figure, red is maximum positive 
energy transfer, while blue is maximum negative energy transfer. Similar plots for r = 830 
(c and d), r = 1000 (e and f), and r = 10000 (g and h). 

FIG. 8: Variation of mean iVTioil and |^^ioi| with e = (r — 1). iVTioil varies as e°'^^. 16*101! 
varies as £"-267 fQj^ e < 1 and as e"'^'^'* for e > 1. 
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FIG. 9: Variation of maximum |n„| as function of reduced Rayleigh number r. 

FIG. 10: Nusselt number vs. e = (r — 1). The Nusselt number shows a power-law scaling 
with e. The power-law exponent is 0.33 ± 0.01. 
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